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ABSTRACT 

Context. Transitional disks are circumstellar disks with dust gaps thought to be related in some cases with planet formation. They 
can shed light on the planet formation process by the analysis of their gas and dust properties. T Cha is a young star surrounded by a 
transitional disk with signatures of planet formation 

Aims. The aim of this work is to study the outer disk around T Cha and to derive its main properties. 

Methods. We have obtained high-resolution and high-sensitivity ALMA observations in the CO(3-2), ^^CO(3-2), and CS(7-6) emis¬ 
sion lines to reveal the spatial distribution of the gaseous disk around the star. In order to study the dust within the disk we have also 
obtained continuum images at 850 pm from the line-free channels. 

Results. We have spatially resolved the outer disk around T Cha. Using the CO(3-2) emission we derive a radius of ~230 AU. We 
also report the detection of the ^^CO(3-2) and the CS(7-8) molecular emissions, which show smaller radii than the CO(3-2) detection. 
The continuum observations at 850 pm allow the spatial resolution of the dusty disk, which shows two emission bumps separated by 
~40 AU, consistent with the presence of a dust gap in the inner regions of the disk, and an outer radius of ~80 AU. Therefore, T Cha 
is surrounded by a compact dusty disk and a larger and more diffuse gaseous disk, as previously observed in other young stars. The 
continuum intensity profiles are different at both sides of the disk suggesting possible dust asymmetries. We derive an inclination of 
i(°) = 67+5, and a position angle of PA(°) = 113+6, for both the gas and dust disks. The comparison of the ALMA data with radiative 
transfer models shows that the gas and dust components can only be simultaneously reproduced when we include a tapered edge 
prescription for the surface density profile. The best model suggests that most of the disk mass is placed within a radius of R < 50 AU. 
Finally, we derive a dynamical mass for the central object of 1.5+0.2 M©, comparable to the one estimated with evolutionary 
models for an age of ~10Myr. 

Key words, stars: pre-main sequence — stars: kinematics and dynamics — stars: individual: T Cha — protoplanetary disks — 
techniques: interferometry 


1. Introduction 


T Chamaeleontis (T Cha) is a young (~7+5Myr) nearby 
(108 pc) T Tauri star in the 6-Cha ass o ciation, surrounded 
by a transition disk (lAlcala et al.l ITo^ iBrown et^ 120071: 
[Torres et ^120081: iMurphv et al .1120131) . There is evidence of a 
dust gap within the disk, a nd a yet unconfirmed substellar com¬ 
panio n inside the gap (see|Huiiamo[et^|201ll;[oloftsonet^ 
l2013l) . If confirmed, the disk around T Cha can give us impor¬ 
tant clues about the physical conditions for substellar formation 
at early evolutionary phases. 


T Cha is surrounde d by a very narrow inn er disk that extends 
from 0.13 to 0.17 AU (lOlofsson et al.ll201^ . and an outer disk 
whose main properties have been inferred from the modelin g 
of its spectral ener gy distribution (SED, e.g. iBrown et~^l2007l) . 
ICieza et ^ (l201ll) showed that the models are highly degenerate 
and can fit the SED of T Cha equally well either with a very 
compact outer dust disk (a few AUs wide) or a much larger but 


more ’tenuous’ disk, with a very steep surface density profile. 
Both family of models suggest a very peculiar outer disk with 
little or no dust beyond ~ 40 AU. On the oth er hand, the cold ga s 
in the T Cha outer disk has been studied bv ISacco et al.l (l2014l) . 
Their spatially unresolved observations suggest the presence of 
a gaseous disk with an outer radius of Rco~80 AU in Keplerian 
rotation. 

Overall, the disk around T Cha shows properties similar to 
the so-called ’faint’ disks, characterized by weak millimeter con¬ 
tinuum emi ssion that can be r esult of different properties or pro¬ 
cesses fe.g. iPietu et al.ll20T4l) . In the case of T Cha there is evi- 
dence of dust clear i ng, grain growth, and a hi gh disk inclination 
(IBrown et al .1120071: iPascucci & Sterzil3l2009l) . 

In this work we present high quality observations of T Cha 
obtained with the Atacama Large Millimeter Array (ALMA), 
which have allowed us to spatially resolve the outer disk around 
T Cha for the first time. ALMA has allowed us to derive basic 
parameters of the outer disk, to break the degeneracy of radiative 
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transfer models based on SED fitting, and to understand if it is 
peculiar in comparison with other circumstellar disks. 


2, Observations 

The observations were performed on 2012 July 01, 26 and 
November 03 at Band 7, as part of the ALMA Cycle 0 program 
2011.0.00921.S. The field of view was ~ 18". A total of three 
data sets were collected, using between 18 and 23 antennas of 
12 m diameter and accounting for 6 hours of total integration 
time including overheads and calibration. Weather conditions 
were good and stable, with an average precipitable water vapor 
of 0.7 mm. The system temperature varied from 150 to 250 K. 

The correlator was set to four spectral windows in 
dual polarization mode, centered at 345.796 GHz (CO(3-2)), 
342.883 GHz (CS(7-6)), 332.505 GHz (802(4(3, l)-3(2,2))), 
and 330.588 GHz (^^CO(3-2)). The effective bandwidth used 
was 468.75MHz, providing a velocity resolution of ~ 0.11 
kms“^ after Hanning smoothing. 

The ALMA calibration includes simultaneous observations 
of the 183 GHz water line with water vapor radiometers, which 
measure the water column in the antenna beam, later used 
to reduce the atmospheric phase noise. Amplitude calibration 
was done using Juno and Titan, and quasars J1256-057 and 
J1147-6753 were used to calibrate the bandpass and the com¬ 
plex gain fluctuations respectively. Data reduction was per¬ 
formed using version 4.1 of the Common Astronomy Software 
Applications package (CASA). We applied self-calibration us¬ 
ing the continuum and we used the task CLEAN for imaging the 
self-calibrated visibilities. The continuum image was produced 
by combining all of the line-free channels using uniform weight¬ 
ing (synthesized beam 0.52" x 0.34", P.A. ~29°; rms = 0.7 mJy 
beam“^). Eor the CO(3-2) line we used Briggs weighting (beam 
0.64" X 0.48", P.A. ~31°; rms per channel = 9 mJy beam“^), and 
for the rest of the lines we used natural weighting, providing a 
synthesized beam -0.8" x 0.6", P.A. -25° and an rms per chan¬ 
nel of 11 mJy beam“^ for the ^^CO and 7 mJy beam“^ for the CS 
and the SO 2 . 


3. Results and discussion 

3.1. Molecular emission line detections 


sion is spatially resolved, with a radius Ruqq - 170 AU after 
deconvolution with the synthesized beam. 

The radial intensity profiles of the CO(3-2) and ^^CO(3-2) 
transitions are displayed in PigureO They have been computed 
using slices along the semi-major axis of the disk (dashed white 
line in Pig.l, left panel) in the gas emission maps. We do not see 
a significant difference of the profiles at both sides of the disk 
(NE and SW) suggesting a symmetric distribution of the gas. 

Gas emission from CS(7-6) is also spatially resolved and 
shows a deconvolved radius of Rqs ~ lOOAU. The integrated 
intensity emission above 3cr is 0.54 + 0.08 Jykms“^. To our 
knowledge, this is the first time that such a high transition of the 
CS molecule is spatially resolved in a disk around a late-type star 
(the first detection in a H erbig star has been recently reported by 
Ivan der Plas et al.ll2014l) . Su lfur-bearing molecules (like CS and 
SO 2 ) studies are scarce fe.g. lDutrev et al.|[201ll:lGuilloteau et al.l 
l2012h . but interesting since these type of mol ecules are present 
in come ts (e.g. Hale-Bopp: [I^da et al.l [20021 and Shoemaker- 
Levv/9: iMatthews et al. 1120021) and thdr presence in the inner 
part of the disks offers the possibility of studying the chem¬ 
ical composition of the planet-forming regions. So far, only a 
few CS detections at lower frequency transi tions have been re¬ 
ported in disks surr ounding K5 to MO stars (iDutrev et al.ll201 lb 
iKastner et al.ll2013h . Here we show for the first time a spatially 
resolved detection of the CS(7-6) transition around a KO star. 

The gas detected in the disk surrounding T Cha shows a 
velocity profile consistent with a Keplerian rotation pattern, 
which can be seen in all the molecular transitions observed (see 
PigureO. CO(3-2) emission is detected at velocities between 
-5.0 to 16.5kms“\ ^^CO(3-2) between -3.0 and 15.0 kms“^ 
and the CS(7-6) between 0.0 to 11.0 kms“^ In Pig|4]we have 
represented a position-velocity diagram along the major axis of 
the gaseous disk traced by CO(3-2). The systemic velocity is 
5.95+0.22 kms“^ In that plot an absorption feature at veloci¬ 
ties near 4.7kms“^ can be seen, which is very likely produced 
by a f oreground molecula r cloud in the same line of sight as 
T Cha (iNehme et al.ll2008l) . By comparing the position-velocity 
diagram with a Keplerian velocity profile (Pigure 11, we show 
that the emission is compatible with Keplerian rotation around 
a 1.3 - 1.7 Mo object in a disk inclined at 67 degrees. This 
mass range is in good agreement with est imations from evolu¬ 
tionary tracks for an age of -lOMyr fe.g. lSchisano et al.]l2009l: 
iMurphv et al.ll2Ql^ . 


Molecular line emission is detected for the transitions CO(3-2), 
^^CO(3-2), and CS(7-6). All of them are spatially resolved for 
the first time for T Cha. No detection was found for 802(4(3,1)- 
3(2,2)), with a 3cr upper limit of 20 mJy. 

Pigure [T] shows the integrated emission maps of the three de¬ 
tected molecules while Pigure |2] displays the intensity-weighted 
mean velocity maps. The CO(3-2) emission is spatially well 
resolved along the major and the minor axis. The integrated 
intensity averaged over the whole emission area above 3cr is 
12 .46 + 0.11 Jy kms ~\ in agreement with the value reported 
by ISacco et ^ (l2014l) from the fit with a Keplerian disk model 
profile to single-dish observations. Considering contours above 
3(T, the outer radius extends to 2.1” which, after deconvolu¬ 
tion with the synthesized beam, corresponds to an outer radius 
Rco ~ 230 AU. The inclination (/) and the position angle (PA) 
of the gaseous disk have been estimated by fitting a Gaussian 
to the CO(3-2) integrated emission map, providing values of 
/ = 67±5°andPA= 113 + 6°. 

^^CO(3-2) emission line shows an integrated intensity of 
4.31+0.07 Jy kms“^ (above 3cr contour). The region of emis¬ 


3.2. Continuum emission at 850 jum 

Continuum emission at 850yum is detected centered at the posi¬ 
tion R. A. (J2000) = 11^57"^13H2,Dec(J2000) =-79°21'31"696. 
The fiux density integrated over the disk structure above 5(T is 
S 850 yum = 198 + 4 mJy. 

Figure \T\ shows the dusty outer disk around T Cha repre¬ 
sented by black contours. The disk is spatially resolved in its 
major axis with a projected diameter of 1.50" (measured at the 
5cr contour level), which corresponds to an outer dust disk radius 
of Pout ~ 80 AU after deconvolution with the synthesized beam, 
and adopting a distance of 108 pc. Two local peaks are observed 
at a projected separation of 0.37” (40 AU at 108 pc), which is 
close to the beam size, and suggest the presence of a gap in the 
inner regions of the disk as predicted by SED modeling. 

In Figure [ 3 ] we have represented the continuum radial inten¬ 
sity profiles at both sides of the disk, together with the average 
profile, including only datapoints above 5-cr. As in the case of 
the gas molecules, they have been computed using slices along 
the semi-major axis of the disk (dashed white line in Pig.l, left 


2 


























Huelamo et al.: ALMA observations of T Cha 


(Jy/beam Mm/ s) 

0,5 1.5 2 2,5 3 


o 

o 


C0(3-2) + continuum 



11^57*^14^0 13®,4 13®.0 


J2000 Right Ascension 


13C0(3-2) + continuum 



11^57*^14®,0 13^4 13®,0 

J2000 Right Ascension 


CS(7-6) + continuum 
30" 



0 0 


1l''57*^14®,0 13®,4 13®,0 

J2000 Right Ascension 


Fig. 1. Integrated emission maps of the CO(3-2), ^^CO(3-2), and the CS(7-6) transitions (from left to right). The black contours 
represent the continuum emission at 850 yum at 5, 15, 30, 45, 60, 75, 90, and 110 cr where 1 cr is 0.7 mJy beam“^ We detect 
two emission bumps separated by 40 AU and an outer dust radius of 79 AU. The white ellipses are the synthesized beams for the 
spectral emission lines and the green ellipse is the synthesized beam for the continuum map. The white dashed line in the left panel 
represents the axis where the position-velocity diagram in Figure |4| has been obtained. 
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Fig. 2. Intensity-weighted mean velocity maps (first-order moment, 2(t cut for CO(3-2) and ^^CO(3-2), and 1.5(T cut for CS(7-6)). 


panel) in the continuum emission maps. We can see a signifi¬ 
cant difference between the profiles at both sides, with the NW 
side being slightly larger than the SE one, suggesting asymme¬ 
tries within the dusty disk. Finally, the i and PA values that we 
derive using the continuum emission at 850 jjm are similar to 
those estimated with the CO(3-2) observations. 


3.3. Comparison with radiative transfer models 


Our ALMA observations reveal that T Cha is surrounded by 
a compact dusty disk with a sharp outer edge at ~80AU and 
a larger gaseous disk with an outer radius of -230 AU. This 
trend, a compact dust disk with a larger and more diffuse 
gaseous disk, has already been observed in a significant num¬ 
ber o f circumstellar disks (e.g. Ilsella et al.|[2007l: iHughes et alJ 
20081: [Andrews et all 120121: Ide Gregorio-Monsalvo et al.l 120131: 
Pietu et al.ll2014l). 

ICieza et all (l201 ll) us ed the radiative transfer code MCFOST 
(iPinte et all 120061 l2009h to model the SED of T Cha. They 
showed that there are two families of dust disk models that can 
reproduce equally well the SED: very small disks (a few AUs 
width) or much larger (/?out~300 AU) but with a very steep sur¬ 


face density profile (with an exponent of a < -2, for a power- 
law prescription). The large degeneracy between these two disk 
parameters. Rout an d a, did not allow the m to choose between 
these two scenarios. lOlofsson et al.l (l2013l) used MCFOST to fit 
the SED together with near-IR interferometric data. They fixed 
Rout=25 AU, and a=-l, and found a best fit model with a narrow 
dust outer disk (Rin=12+2 AU). 

Our ALMA data shows that models with a power-law surface 
density, like the ones used in these two works, cannot reproduce 
both the CO and continuum data. We nevertheless first consider 
this type of models to discuss the ALMA data in the context of 
previous results. 

We have modeled the SED of T Cha with M CFOST, be¬ 
ing th e starting point the model grid prese nted bv ICieza et al.l 
(1201 ll) and refined bv lOlofsson et al.l (l2013l) . Basically, the disk 
is composed of 2 sub-disks: an inner and an outer disk with a 
density structure defined by a power law surface density pro¬ 
file with exponent a, S(r) = Eq (r/ro)^, and a scale height of 
h(r) = with (3 being the disk flaring index, and /to 

the scale height at a reference radius vq = 50 AU. Each disk ex¬ 
tends from an inner radius, Riu to an outer radius, Rout- The grain 
size distribution in each disk is defined by &n{a) oc a^da, be- 
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Fig. 3. Top and middle panels: Continuum and gas radial inten¬ 
sity profiles. The blue (SE) and red (NW) data points are the 
ALMA data (over 5 -(t) at both sides of the disk. In the case of 
the continuum we have also included the average profile (green 
data). The solid black lines show the best model using a tapered 
edge prescription for the surface density, while the dashed lines 
show the best model using a truncated power law. The upper 
panels show the continuum data, including a zoom of the outer 
regions, while the middle panels include the CO profiles. Bottom 
panel: the observed SED and the fit from the tapered edge model 
presented here. 


tween a^nm and timax- The temperature structure in the disk is 
calculated by considering the dust opacity only. The dust prop¬ 
erties are computed assuming the Mie theory. Einally, the total 
dust mass in the disk (including all grain sizes) is represented by 

^dust- 

To calculate the CO channel maps and surface brightness dis¬ 
tribution, we assume a constant gas-to-dust mass ratio of 100 
throughout the disk (both radially and vertically). We adopted 
a standard CO abundance with respect to H 2 (10“"^), set con¬ 
stant through the disk where rdust> 20 K and equal to zero where 
7'dust< 20 K to mimic the effect of CO freeze out. The ^^CO/^^CO 
ratio is set to 76. The level populations are calculated assuming 
T'gas = T'dust at each point in the disk . The radial and vertical tem¬ 
perature profiles and the radiation field estimated by the Monte 
Carlo simulation are used to calculate level populations for the 
CO molecule and to produce the SED, continuum images, and 
line emission surface brightness profiles, as well as kinematics 
with a ray-tracing method. The kinematics are calculated assum¬ 
ing the disk is in pure Keplerian rotation. 

We have adopte d the inner disk parameters from 
lOlofsson et all (l2013l) . Eor the outer disk, we have fixed 
two parameters obtained from the ALMA data, /?out and /, se¬ 
lecting the grid values closer to the ALMA measurments (80 AU 
and 68°, respectively). We have explored /?in, P, /to@50AU, 
tr, and Mdust, using the same parameter range shown 
in Cieza et aO (l201 ll) . Eor the SED, we have fitted the same 


CO(3-2) 



- 2-1012 
Position offset (arcsec) 


Fig. 4. Position-velocity diagram along the disk major axis (see 
exact axis represented in Eig[T]left panel). The dotted, solid, and 
dashed black curves represent the best fits to the data, that cor¬ 
respond to a Keplerian velocity profile for a systemic velocity 
of 5.95 km s“\ an inclination of /(°) = 67, and a central mass of 
1.3, 1.5, and 1.7 M©, respectively. 


observational dataset displayed in that work. The ad opted stellar 
parameters are reff=5400K , Av=1.5, and J=108pc (iTorres et al.l 
l2008l:ISchisano et al.ll2009l) . The best disk model, that is, the one 
with the minimumprovides parameters of a= -2.5, 1.07, 

/zo@50AU = 6 AU, /?in = 19 AU, = 1000, and a disk dust 
mass of Mdust -1x10“^ M©. 

Since SED modeling is highly degenerate, the best-fit model 
is unlikely to be a unique solution. Therefore, we have performed 
a Bayesian analysis to estimate the validity range for each of the 
explored parameters (iPress et al.lll992l: iPinte et aP l2007l) . The 
result is displayed in Eigure [3 where we show the Bayesian 
probability distributions for the different disk parameters. While 
Afdust, ^in and ho seem well constrained, this is not the case for a\ 
it shows a local peak at a = -2.5 but a relatively fiat distribution. 
We conclude that, even fixing /?out, remains unconstrained by 
the SED modeling. 

With our ALMA observations we have partially broken the 
(of, /?out) degeneracy commonly encountered with SED fitting by 
measuring /?out- The resolution reached by our observations does 
not allow us to constrain accurately the surface density profile. 
But with Rin ~ 20 AU and a disk width of ~60 AU, we exclude 
surface density profile shallower than -1. According to the ob¬ 
served Rout, we can also discard the family of models with very 
narrow dusty rings, and the extreme case of a very large disk 
(Rout ~300 AU) with a —3. 

If we take the gas emission into account, the model fails to 
fit simultaneously the gas and dust profiles (see Eigure O, as 
alre ady observed in othe r spatially resolved circumstellar disks 
fe.g. iHughes et al.ir2008h . As discussed by the authors, a power 
law density profile cannot reproduce the different extent of the 
gas and dust emission observed in circumstellar disks while 
a tapered edge model, in which the surface density falls off 
gradually, can in principle reconcile the observed profiles. We 
have therefore run the MCEOST code, but using a tapered ex¬ 
ponential edge in the surface density distribution. In this case, 
the density profile is represented by a function defined by the 
characteristic radius, Rc (the radius out of which the brightness 
drops towards zero) and 7 , the surface density index: S(r) = 
SoWronexp(-(r/ro)2A. 
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Fig. 5. Bayesian probability distributions of the disk parameters 
of T Cha. The blue lines show the results without fixing any 
parameter, while the red lines show the results after fixing the 
outer disk radius and the disk inclination. 


For t his model, we have als o adopted the inner disk parame¬ 
ters from lOlofsson et al] (l2013l) . and varied only the prescription 
for the outer disk. We have fixed the disk inclination to 68°, the 
dust mass to Mdust= 9x 10“^ M© (de Gregorio-Monsalvo et al., in 
prep.), and we have explored the rest of the disk parameters: we 
have sampled a range of y values between 0 and 2, Rc between 
40 and 100 AU, Rin between 12 and 30 AU, the scale height at 
50 AU (/zo @50AU) between 3 and 6 AU, and the fiaring index, yS, 
between 1.0 and 1.1. Finally, to take into account the difference 
in brightness between both sides of the disk in the continuum, 
we explore models that pass between the two observed profiles. 

The result is displayed in Figure [3] where we show the best 
model that can fit simultaneously the two disk components and 
the observed SED. The model shows y = 0.5 and Rc = 50 AU for 
both the gas and the dust. We also derive these parameters: Rin= 
20 AU, /zo@50 AU = 4 AU, and p = 1.0. This model is consis¬ 
tent with having the gas and the dust well mixed an d mainly lo¬ 
cated a t a radius smaller than 50 AU, as suggested bv ICieza et al.l 
(I 2 OIII) based on the steep drop of the SED at sub-mm wave¬ 
lengths. 

Figure [ 3 ] shows that our best model does not perfectly fit the 
CO line profiles, which can be related either with the underly¬ 
ing chemistry (we assumed ISM abundances and very simple 
CO freeze-out) or with the model prescriptions. In fact, tapered- 
edge models sometimes fail to reproduce simultaneously the ob¬ 
served dust and gas profiles obtained from very high spatial res¬ 
olution and sensitivity observatio ns (see [Andrews et al.l l2Q12l: 
Ide Gregorio-Monsalvo et al.ll2Q13l) . and suggests that other pro¬ 
cesses like e.g. grain growth and radial migration should be taken 
into account. Given that the disk is barely resolved in our obser¬ 
vations, we expect future, higher spatial resolution observations 
to provide stronger constraints on the relative location of the gas 
and dust, and the departure from point-symmetry. 

4. Conclusions 

High spatial resolution and high sensitivity ALMA observations 
have allowed us to spatially resolve the outer disk around the 
young and isolated object T Cha. The target is surrounded by 
a compact dusty disk and a ~3 times larger gaseous disk. Our 
main results can be summarized as follows: 

- We have spatially resolved the gaseous disk of T Cha in three 

different molecular emission lines: CO(3-2), ^^CO(3-2) and 


CS(7-8). Using the CO(3-2) image we derive an outer ra¬ 
dius of /?gas,out -230 AU, an inclination of /(°) = 67+5, and a 
position angle of FA(°) = 113+6. The line intensity profiles 
are similar at both sides of the disk in the CO molecules, 
consistent with a uniform distribution of the gas. 

- The disk around T Cha is in Keplerian rotation, and 
the estimated dynamical mass of the central object, 

= I. 5 +O. 2 M 0 , is in good agreement with previous esti¬ 
mations based on evolutionary tracks. 

- The dusty disk is resolved in the continuum observations 
at 850yum and it shows a similar i and PA to the gaseous 
disk. The continuum intensity profile displays two emission 
bumps separated by 40 AU, suggesting the presence of an 
inner dust gap as predicted by SED modeling, and an outer 
radius of ~80 AU. The profiles are different at both sides of 
the disk, which points towards asymmetries in the dust dis¬ 
tribution. These data allows us to rule out both the very small 
and large /?out families of SED models. 

- Radiative transfer models including a truncated power law 
prescription for the surface density profile cannot reproduce 
simultaneously the gas and dust profiles. We can fit both 
components simultaneously using a tapered-edge model pre¬ 
scription for the surface density. The best model provides 
values of y= 0.5, and Rc= 50 AU, which is consistent with 
having most of the disk mass within the inner 50 AU. 
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